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ABSTRACT 


— ■>  The  difference  between  consecutive  order  statistics  from 
a  sample  is  called  a  spacing.  Various  tests  based  on  sample 
spacings  have  been  considered  in  the  literature  for  testing 
the  hypothesis  that  the  sample  is  drawn  from  a  specified  dis¬ 
tribution.  Tests  based  on  the  spacings  arc  recommended  for  use 
when  the  alternative  distribution  differs  from  the  hypothetical 
distribution  in  the  shape  of  the  density  function.  In  this 
paper,  we  consider  a  test  based  on  the  spacings  designed  for 
the  case  when  the  ratio  of  the  two  density  functions  i 
piece-wise  monotone  function.  This  paper  deals  main 
the  large  sample  properties  of  the  test.^^^^^ 
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1. 


Introduction. 


Let  X^,...,X  be  a  sample  drawn  from  a  cer¬ 
tain  distribution.  In  this  paper,  we  consider  a  Lost  of  the 
hypothesis  Hq  that  the  sample  comes  from  a  known  distribution 
F,  say,  against  the  alternative  hypothesis  1!^  that  the  sample 
is  drawn  from  a  distribution  G,  say,  where  G  is  not  known  com¬ 
pletely.  Let  A  denote  the  common  support  of  the  distributions 
F  and  G.  We  assume  that  A  is  a  finite  or  infinite  interval  and 
that  distributions  have  continuous  density  with  respect  to  the 
Lebesgue  measure.  Let  f  and  g  denote  the  respective  density 
functions.  We  further  assume  that  the  ratio  Q(x)  =  g(x)/f(x) 
is  a  piece-wise  strictly  monotone  function  of  x  inside  A. 

Suppose  that  the  slope  of  the  graph  of  0(x)  changes  sign  at 
k  points.  Let  n  =  ( , . . . , n^)  denote  the  change  points.  We 
consider  two  cases:  (i)  _n  is  known  and  (ii)  jj  is  not  known. 

We  shall  consider  in  length  the  case  in  which  there  is  a 

single  change  of  sign  (k=l)  in  the  slope  of  the  graph  of  Q(x). 

This  is  realized  in  the  following  examples,  where  <Mx)  = 

-1/2  2 

(2"’)  exp(-x  / 2)  denotes  the  standard  normal  density  function. 


(a) 

f  (x) 

=  ■'  (x)  ; 

g(x)  = 

(l/°)  V  (x/o) 

(b) 

f  (x) 

=  4>  (x)  ; 

g  (x)  = 

e"x/ (l+e_X) 

(c) 

f  (x) 

=  Mx); 

g  (x)  = 

Pv  ( X-  u  )  +  (  1 

VI1 

<N 

V 

o 

(d) 

f  (x) 

=  1,  O'x'l; 

g  (x)  = 

I  ('(  +  !')  if 

](<>)!  (1  ) 

0  <  x  -■  1; 


n  ,  ('■  1  . 
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Lot  =  F(X^),  i  -  1 , . . .  , n  and  lot  D ^  -  Y  ^  , 

D2  =  Y  (2) -Y  ( i)  '  *  *  *  ,Dn  =  Y(n)_Y(n-l)  dcnott!  thc  sample  spacings 
of  the  transformed  data  whore  Y^  denotes  the  ith  smallest 
value  among  Y  ,...,Yn.  Our  tost  is  based  on  the  spacings. 

The  test  statistic  will  be  denoted  by  T.  Various  tests  for 
goodness-of- f it ,  based  on  spacings  have  been  considered  in 
the  literature.  The  papers  of  Pykc  (1965,  1972)  ,  Proschan  and 
Pyke  (1967)  ,  Sethuraman  and  Rao  (1970)  ,  Kale  (1969)  and 
Kirmani  and  Alam  (1974)  may  be  cited  for  reference.  Pyke  (1972) 
points  out  that  tests  based  on  spacings  should  be  used  when  the 
alternative  distribution  differs  from  the  hpothetical  distri¬ 
bution  in  the  shape  of  the  density  function. 

The  following  statistics  have  been  proposed  in  the  litera¬ 
ture  for  a  test  of  goodness-of- fi t: 


U 


n 

r 

i=  1 


1 

2  ' 


n 

V  =  >  | nD . -1 |  and 

i  =  l  1 

n 

W  =  )  log  D . . 

i  =  l 


Then  null  hypothesis  is  rejected  when  thc  absolute  value  of  the 
statistic  is  large.  It  is  known  (see  e.g.,  Cibisov  (1961)) 
that  the  asymptotic  efficiency  of  any  test  symmetric  in  the 
spacings  is  equal  to  zero  relative  to  the  Kolmogorov-Smirnov 
test.  Sethuraman  and  Rao  (1970)  have  compared  the  relative 
efficiencies  of  the  tests  based  on  U,  V  and  W. 
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For  applications  to  reliability  and  li fe-testing ,  Proschan 
and  Pyke  (1967)  have  considered  a  test  of  the  hypothesis  that 
the  given  sample  comes  from  an  exponential  distribution  which 
has  a  constant  failure  rate  property  against  the  alternative 
hypothesis  that  the  distribution  has  monotone  failure  rate. 

The  test  is  based  on  the  statistic 

n  n 

(1.1)  S  =  l  l  h (D .  , D  . ) 

i  =  l  j  =i  1  J 

whore  h  is  a  bounded  nonnegative  function  and  =  (n-i+l)D. 
denotes  the  ith  normalized  spacing.  The  authors  have  shown 
that  the  distribution  of  S  is  asymptotically  normal  under  the 
alternative  hypothesis.  Further,  Bickel  and  Doksum  (1969) 
have  shown  that  the  asymptotic  normality  holds  also  for  a 
sequence  of  alternatives  approaching  the  exponential 
distribution . 

The  test  statistic  T  is  derived  from  S  as  follows: 

First  suppose  that  r  is  known.  Let  =  F(n.)  and  let  ~  . 

-  ii  i 

denote  the  set  of  values  of  Y.  for  which  .  ..  ■'  Y.  •  ■  .  ,  and 

3  i-l-3i 

n.  denote  the  number  of  elements  in  n . ,  i  = 
i  1  ' 

* 

where  '  =  0  and  i,  ,  =  1.  Let  Y.  denote  the  smallest  value 

o  k+1  i 

* 

in  .  ,  and  let  S  =  Y .  -  . .  ..  for  m  =  n.  +  ...+n.  .+1  (i  =  l,...,k+l) 

l  mi  l-l  1  l-l 

and  S  =  D  otherwise,  m  =  l,....n,  where  n  =0.  Let 
mm  o 

1  for  x  £  y 

(1.2)  h (x , y )  =  { 

0  for  x  >  y . 


Define 
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T  T  or  T 
o  o 

according  as  the  slope  is  positive  or  negative.  The  hypothesis 

H  is  rejected  for  small  values  of  T. 
o 

In  case  (ii)  ,  we  estimate  _n  from  the  data  and  substitute 
the  estimate  for  g  in  the  definition  of  T,  given  above.  The 
estimation  of  is  considered  in  Section  5  below. 

In  Section  2,  we  show  that  the  given  test  is  unbiased. 

To  compute  the  critical  value  of  T  and  the  power  of  the  test, 
we  need  to  find  the  distribution  of  T  under  Hq  and  H^.  The  dis 
tribution  is  shown  in  Sections  3  and  4.  Some  results  on  the 
relative  efficiency  of  the  test  are  given  in  Section  6. 
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2.  Unbiasedness  of  the  test.  Let  11  =  GF  ^ =  H(Y^)  and 

D  !  =  Y ,  -  Y .  ,  ,  i  =  1 , ,n  ,  where  Y '  =  0  .  Suppose  that 
1  (i)  (l-l)  o 

the  slope  of  the  qraph  of  Q(x)  is  positive  inside  (<i  ,  ,  ••  ). 

in- 1  m 

Then  H  is  convex  on  (f  ,  ,>~  ).  Therefore,  by  the  mean  value 

m- 1  m 

theorem  we  have  that  for  n,  +  ...+n  ,  +  l  •  i  <  i  ■  n,  +  ...+n 

1  m- 1  —  —  -  1  m 


-  Y ' 

(j-1) 


-Y\ 


(i-1) 


v  -  Y 

-ill  (j-1) 

Y  -Y 

<i>  (i-1) 


H’  (C.) 
H'  (Ci) 


(2.1) 


Y  -Y 

-ill  (i-1) 

Y  -Y 

(i)  (i-1) 


D  . 

D. 

l 


where  Y,.  ..  ■  C .  Y,..  and  Y,.  . .  ■  C.  <  Y,...  The  inequality 

(]-l)  3  (:)  (i-l)  i  (i) 

in  (2.1)  follows  from  the  fact  that  H'  (C ^ ) /H *  (CJ  _  1  due  to 

the  convexity  of  11  on  (t  .  , r  ).  From  (2.1)  we  find  that 

D.  '  D.=>  D'  <  D'..  Therefore  h(D.,D.)  <  h(D'.  ,D'.).  Let  T’,  T’ 
1-3  1-3  13-13  mm 

and  T'  be  obtained  from  T  ,  Tm  and  T  respectively  by  substituting 

Y f  for  Y^  ,  i  =  l,...,n.  Then  by  the  above  inequality  we  have 

that  T  '  T' .  Similarly,  if  the  slope  of  the  graph  of  0 (x) 
m  —  m 

is  neqative  inside  (n  - , n  )  then  H  is  concave  on  (f  ,  ) 

^  m-1  m  m-1  m 

and  therefore  T  <  T' •  It  follows  that  T  <  T' .  Hence 
m  —  m  — 


P  (T  _  t !  G)  _  P (T 1  £  1 1 G)  =  P (T  £  t|F),  for  all  t  _  0. 


Since  H  is  rejected  for  small  values  of  T,  the  test  is 
unbiased . 

3.  Distribution  of  T  under  H  .  When  G  =  F,  the  sample  spacings 
'IV' 

t  , . . . ,Dn  are  jointly  symmetrically  distributed  according  to 

the  Dirichlet  distribution  given  by  the  density  function 


1 


iur 


(3.1)  p(d.,...,d  )  =  n!  ;  d  '•  0,  0  ■'  d.+...+d  '  1. 

1  n  l  ~  —  1  n  — 

Let  Z^,...,Z  be  n  random  variables  jointly  symmetrically 
distributed  and  let 


R  . 


i=l 


h (Z  i  ' Z  j  )  '  3  =  1 /  ■  •  •  /  n 


denote  the  left  sequential  ranks  where  the  function  h  is  defined 
as  in  (1.2).  Similarly,  let 


n 

V 


R.  =  i  h(Z.  ,2  .)  ,  j  =  1,.  .  .  ,n 
3  i= j  13 


denote  the  right  sequential  ranks.  It  is  known  (see,  e.g., 
Renyi  (1962))  that  R^,...,R  are  statistically  independent  and 
that  the  distribution  of  is  uniform,  given  by 


(3.2) 


P(R.  =  m)  =  -r 


m  —  l,...,j. 


The  same  result  holds  for  R.,...,R  . 

1  n 

The  above  result  is  easily  generalized  as  follows:  Let 

1  1  1  +  1  1+  2  V  +V1 . l+’-*  +  \+l 

be  a  partition  of  Z  =  (Z^,...,Z  )  into  k+1  sub-vectors,  where 

■  ,  +  .  =  n  and  let 

1  k+1 


=  ?  h(2,  „  ... 

i  .  u  t  ,  +  •  •  •  r  it. 


i=l 


T+i'2 ?.  +  ■••  + p  ,  +  j  ^ 
m-1  1  m-1  J 


be  the  jth  left  sequential  rank  of  the  variables  in  the  mth 
partition  (Z 


]?.,  +  •••+?.  .  + 1 
1  m-1 


Zi  m);  3 

1  m 
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m  =  l,...,k+l,  where  =  0.  Let  the  right  sequential  ranks 

R™  be  defined  in  analogous  manner.  Then  the  random  variables 

R™,  j  =  l,...,'m,  m  =  l,...,k+l  are  statistically  independent’ 

and  the  distribution  of  each  R™  is  uniform.  Similarly,  the 

random  variables  R.,  j  =  l,...,1  ;  m=  l,...,k+l  are  statis- 

3  J  m 

tically  independent  and  the  distribution  of  each  R™  is  uniform. 
Moreover,  R™  and  R™  are  independent  for  m  4-  m'  .  These  results 
follow  from  the  property  of  symmetry  of  the  joint  distribution  of 

Zl"  •  *  ,Zn‘ 

Given  £  and  n  =  (n ^  ,  . . .  ,n^+ ^ )  ,  ■'  .  represents  a  sample  of 

n.  observations  from  a  uniform  distribution  on  (f,. 
l  l- 1  l 

Moreover,  the  k+1  subsamples  are  conditionally  independent. 

The  following  properties  of  the  conditional  distribution  of 
T  ,...,T^+^;  T^  , . . .  ,T^+^ ,  given  n  follow  from  the  results  given 
above:  (a)  T^,...,T^  +  j  are  independent,  (b)  Tj  , . . . ,T^+ ^  are 

independent,  (c)  T^  and  T ^  are  independent  for  i  ^  j  and 
(d)  T^  and  T^  have  the  same  distribution  for  each  i.  Therefore 
conditionally , 


(3.3) 


.  d  T. 

T  -  '  x 


where  -  means  "distributed  as". 

Lot  P.  (t)  =  P  (T .  =  t  j  n )  ,  i  =  1 , . . . , k  +  1 .  Clearly, 

' 11  i 

P.  (t)  -  0  for  t  n.  .  For  t  '■  n.  ,  the  probability  can  be 
l  ,n .  i  x  1 

l 

computed  recursively  from  the  relation 


(3.4)  "  P  Ct>  *  L 

1 


^"i-l 


(t-D+P, 


(t-2) +. . . +P. 


(t-n.) 

x,ni— 1  x 
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Since  T  ,  .  . .  ,T  are  conditional ly  independent  given  n, 

J.  K  '  J. 

from  (3.3)  we  get 

k+1 

(3.5)  P(T  =  T|n)  =  l  "■  P.  (t  )  . 

V-+tk+l  1  =  1 

As  n  is  distributed  according  to  the  multinomial  distribution 
with  the  associated  probability  vector  p  =  (p  ,  ...,p,  ), 

p^  =  ^ ,  i  =  l,i..,k+l,  we  obtain 

n  . 

k+1  p.  1 

P(T  =  t)  =  ;  P  (T  =  tin)  n!  ii  — — r  . 

n1+...+n^+1  -  n  i=l  ni ‘ 

Kendall  (1938)  and  Mann  (1945)  have  tabulated  the  distribution 

ofp.  (t),forn.  '-'10. 
l 

Since  Tt  is  the  sum  of  independent  random  variables 

Rw-.-.U1',  it  follows  from  Liaponouv's  theorem  that  for  large 
In. 

x 

values  of  n.,  T ^  is  asymptotically  normally  distributed  with 

moan  n^(n^  +  3)/4  and  variance  n^(n^-l)  (2n^+5)/72.  Therefore, 

T  is  asymptotically  normally  distributed  for  large  n,  with 

2 

mean  .  and  variance  a  /  given  by 
n  n 


k+1  k+1  _ 

(3.6)  =  );  n  (n  +3)/4  -  (n  /4)  \  p* 

°  i=l  i=l  1 


(3.7) 


7  k+1  ^  k+1  , 

■0  =  V  n.  (n.-l)  (2n  +5)  (n  /36)  \  p]?. 

i~l  i  =  l 


where  ~  means  "asymptotically  equivalent  to". 
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— 


4.  Distribution  of  T  under  H  .  The  small  sample  theory  of  the 
distribution  of  T  under  is  mathematically  intractable. 
Therefore,  we  consider  the  asymptotic  theory.  Let  Z^,...,Z 
be  a  sample  from  the  exponential  distribution  whose  cdf  is 
given  by  A(z)  =  l-exp(-z),  and  let  Z^.^  denote  the  ith 
order  statistic  in  the  sample.  Similarly,  let  denote 

the  ith  order  statistic  in  a  sample  of  n  observations  from 
the  uniform  distribution  in  (0,1).  It  is  known  that 


(i) 


V  Z  ./  (n-j  +  1)  . 

5*1  : 


Let 


(U)  --  (1-u)  /H'  (H  (u)  )  , 


where  H  =  GF  ^ .  We  have 


(4.1)  nDi  d  n  ( H  I/\(Z(i))-H  1A(Z(i_1))) 

=  n(Z  (i)-z  (i_d  )e_Z/g(G-1  (A  (/.)))  , 


by  moan  value  theorem 


=  n  (Z  (i)  "Z  (i-1)  )  "  (U> 


(n/  (n-i  +  1 )  Zi  \  bm] 


wher<. 


(i-1) 


Z(i)  ,  u  =  A(z)  and  U(._:)  1  c.  <  U(i). 


Though  (  Z  ^  ,  Z  _.  )  and  (m,m)  are  statistically  dependent, 
Proschar.  and  Pyke  (1967)  have  shown  that  the  degree  of  depen¬ 
dence  is  negligible  in  relation  to  the  distribution  of  T  for 


10 


for  large  n.  Under  the  condition  of  independence,  we  have 


Eh  (Di  , D_.  ) 


P  (D 


D  . ) 
1 


■  (  •  )  !>.(■) 

_ i-  L _ 

n-j+1  !n-i+l 


(  1  .) 


-1 


n-j  +  1 


Then  for  large  n ,  we  have 

ET  .  =  E  ;  h  ( D  .  ,  D  . ) 

1  n._1+l-i£jin.  1  1 

(4.2)  n~  i  ]  IA--LLL ll'JJ _  dudv 

:  .  J  H'  ( H~  1  (u)  )  +11 '  (H_1  (v)  ) 

i-  1 


=  n  ...  ,  say. 


From  Theorem  4.2  of  Proschan  and  Pyke  (1967)  it  follows  that 
T.  is  asymptotically  normally  distributed.  The  mean  of  the  asympto¬ 
tic  distribution  is  given  by  (4.2).  The  variance  of  the  asymptotic 

3  2 

distribution  equal  to  n  say,  can  be  obtained  from  formulas 

(4.59)  and  (4.60)  of  the  paper.  We  do  not  give  the  expression  for 
the  variance,  since  it  is  involved.  The  distribution  of  T  is 

2  k+1 

asymptotically  normal  with  mean  equal  to  n  'i  u  ,  and  variance 

equal  to  n  >  a  ■  . 

i  =  l 

5.  Estimation  of  H.  We  have  considered  the  case  when  or 
equivalently  y  is  known,  giving  the  points  of  inflection  in  the 


A 


11 


graph  of  Q(x).  If  _  is  unknown,  we  estimate  i  as  follows.  We 
shall  describe  the  method  of  estimation  when  there  is  a  single 
point  of  inflection  \  ^ ,  though  the  method  carries  through  to 
the  case  in  which  there  are  several  points  of  inflection.  Let 


L 

m 


mi  n  n 

l  >  h(D  .  ,D  )  +  l  V  h(D  D.) 

j=l  i=l  ^  j=m+l  i=j+l  J 


* 

and  let  m  denote  the  value  of  m,  maximizing  (minimizing)  as 
m  varies  from  0  to  1,  if  the  slope  of  Q(x)  changes  sign  from 
negative  to  positive  (positive  to  negative).  If  the  sign  changes 
from  negative'  to  positive  then  r(u)/(l-u)  is  first  increasing 
then  decreasing  as  u  varies  from  0  to  1.  From  the  representation 

(4.1)  of  the  sample  spacings,  it  is  seen  that  the  values  of 

tend  to  increase  then  decrease  as  i  varies  from  1  to  n  and 

* 

therefore  the  value  of  m  maximizing  L  is  approximately  given 

by 

(5.1)  m  --  (  ;  j  <  n  ] 


where  [x]  denotes  the  nearest  integer  value  of  x.  If  the  slope 

★ 

of  Q  (x)  changes  sign  from  positive  to  negative  then  m  minimizes 
.  The  estimate  of  •'  ^  is  given  by  the  largest  value  of  ' 


satisfying  the  relation  (5.1). 


*  -1/2, 


It  can  be  shown  that  m  /n  =  ?'  +  0  (m  ).  Therefore, 

i  P 

the  asymptotic  theory  developed  in  the  previous  section  remains 
valid  when  ^  is  replaced  by  its  estimated  value  in  the  definition 


of  T. 
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6.  Asymptotic  relative  efficiency.  Tn  this  notion,  we  compare 
the  test  based  on  T  with  a  likelihood  ratio  test,  using  the 
criterion  of  asymptotic  relative  efficiency  (ARE)  for  the  com¬ 
parison.  We  consider  below  two  examples  for  the  comparison.  In 
Example  1  we  test  an  exponential  distribution  agains  a  Weibull 
distribution.  In  Example  2  we  test  a  uniform  distribution 
against  a  beta  distribution.  For  a  specified  set  of  alternatives 
indexed  by  ,  say,  the  formula  for  the  ARE  of  a  sequence  of 
tests  (based  on  a  sequence  of  asymptotically  normal  test  statistics 
! i )  against  a  sequence  of  tests  (based  on  the  asymptotically 
normal  test  statistics  it  ))  is  given  by  the  formula  (sec  e.g., 
Gibbons  (1971)) 


(6.1) 


ARE 


lira  i  n  °  , 

n 


u'  (•■'  )  2 
i  o 
n 


n 


where  denotes  the  null  hypothesis,  ii  (v)  and  a_  (u)  denote 
O  T  T 

n  n 

the  limiting  mean  and  variance,  respectively,  of  1  T^ }  ,  n^,  (") 

n  n 

denotes  the  derivative  of  uT  (0)  with  respect  to  0.  The  parametric 

n 

2 

functions  .!  .  (),  y  ( t> )  and  a'  (t)  are  defined  similarly,  as  for 

'  n  n  1  n 


T 


— 

Example  1.  Let  f(x)  =  e  ,  x 
x  0,  f-  •>  1.  it  is  seen  that  Q(x) 
in  x  for  x  £  (>)  for  all  0  >  1.  We 
The  spacing-test  statistic  is 


f)  —  2.  —  x 

-■  0  and  g(x)  =  Ox  e  , 
is  increasing  (decreasing) 
have  r.j_  =  F  ( 1)  =  1-e  1 . 


T 


T1  + 


.  Ui 
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From  (4.2)  we  have,  after  simplification, 


(6.2) 


14m  :  n  (I x+i 2) 


where  1^  and  are  given  by 


(6.3) 


X1  = 


’2z2e  y~2 [ ( 1+y ) log  y  -  (l+z)log  z]  , 

.  2  2,2  aya 

(y  +z  ) 


0<y<z<l 


,  approximately. 


(6.4) 


T2  = 


y2z2e  ^y+Z^ ^yZ [ ( (1+z) /z) log  z- ( ( 1+y) /y) log  y] 

.  2  2.2 
(y  +z  ) 


0<y-:z<  1 


,  approximately. 


The  approximate  values  of  1^  and  X2  given  above  are  obtained  by 
numerical  integration.  From  (3.7)  we  have 


(6.5) 


0,^(1) 


¥<1'2f'l(l-fl))- 


Proschan  and  Pyke  (1967)  have  shown  that  the  likelihood 
ratio  test  rejects  the  null  hypothesis  when  Tw  is  large  where 


n 

(6.6)  T  =  l  (1-X  )log  X. 

W  .  L ,  l  i 

1  =  1 


and 


2  2 

u'  (1)  =  n[(y-l)  +tt  /6],  y  =  .5772,  approximately. 


dydz 


(6.7) 
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(6.8)  ol  (l)  =  nf  ( Y  —  1 2 + 7r  2  /  6  ]  . 

W 

From  (6.1)  the  ARE  is  given  by 

ARE  =  36(I1  +  I2)2/(l-2f.1  (1- f  ))  (  ( V  - 1 )  2  +  -  2  /  6 ) 


=  approximately. 

Example  2.  Let  f(x)  =  1,  0  x  *  1  and  cj(x)  =  'x  '  1 ,  1. 

It  is  seen  that  Q(x)  is  increasing  in  x  inside  the  interval  (0,1). 

Therefore  T  =  T  .  From  (4.2)  and  (3.7)  we  get  after  some 
,  o 

I 

simplification 

2 

;ii  (1)  "  "  •  •  •  1  :  -(.031)n2 

io  ^  2  5  7 

a2  (1)  ’  n3/ 36 . 

o 


The  likelihood  ratio  test  rejects  the  null  nypothesis  when 
n 

T  =  >  log  x.  is  large.  By  direct  computation 

B  i  =  l  x 

’i '  (1)  =  -n  and  u2  (1)  =  n. 

B  B 


From  (6.1)  the  value  of  the  ARE  turns  out  to  be  equal  to  .0345 
approximately . 
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